Direct measurement of general quantum states using weak measurement 
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Recent work [J.S. Lundeen et al. Nature, 474, 188 (2011)] directly measured the wavefunction by 
weakly measuring a variable followed by a normal (i.e. 'strong') measurement of the complementary 
variable. We generalize this method to mixed states by considering the weak measurement of various 
products of these observables, thereby providing the density matrix an operational definition in terms 
of a procedure for its direct measurement. The method only requires measurements in two bases 
and can be performed 'in situ', determining the quantum state without destroying it. 
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The wavefunction ^ is at the heart of quantum me- 
chanics, yet its nature has been debated since its incep- 
tion. It is typically relegated to being a calculational 
device for predicting measurement outcomes. Recently, 
Lundeen et al. proposed a simple and general opera- 
tional definition of the wavefunction based on a method 
for its direct measurement: "it is the average result of 
a weak measurement of a variable followed by a strong 
measurement of the complementary variable [l|, By 
'direct' it is meant that a value proportional to the wave- 
function appears straight on the measurement apparatus 
itself without further complicated calculations or fitting. 
The 'wavefunction' referred to here is a special case of a 
general quantum state, known as a 'pure state.' The gen- 
eral case is represented by the density operator p, which 
can describe both pure and 'mixed' states. The latter in- 
corporates both the effects of classical randomness (e.g., 
noise) and entanglement with other systems (e.g., deco- 
herence) . The density operator plays an important role in 
quantum statistics, quantum information, and the study 
of decoherence. Because of its generality and because 
it follows naturally from classical concepts of probabil- 
ity and measures, some consider p, rather than "5 , to be 
the fundamental quantum state description. In this let- 
ter, we propose two methods to directly measure general 
quantum states, one of which directly gives the matrix 
elements of p. 

The standard method for experimentally determining 
the density operator is Quantum State Tomography Q. 
In it, one makes a diverse set of measurements on an 
ensemble of identical systems and then determines the 
quantum state that is most compatible with the mea- 
surement results. An alternative is our direct measure- 
ment method, which may have advantages over tomog- 
raphy, such as simplicity, versatility, and directness. A 
quantitative comparison of measures such as the signal 
to noise ratio, resolution, and fidelity, has not been un- 
dertaken but some limitations of the direct method have 
been identified in Q . As compared to tomography, which 
works with mixed states, the most significant limitation 
of the direct measurement of the wavefunction is that it 
has only been shown to work with pure states. 



Previous works have developed direct methods to mea- 
sure quasHjrobability distributions, such as the Wigner 
function [5j, Husimi Q- function [6[, and the Glauber- 
Sudarshan P-function [7| • These are position-momentum 
(i.e. 'phase-space') distributions that are equivalent to 
the density operator, and have many, but not all, of the 
properties of a standard probability distribution. The 
Wigner function can be directly measured by displacing 
the system in phase space and then measuring the parity 
operator Q. Equivalently, the integral of the interfer- 
ence between a pair of rotated and displaced replicas of 
the system will give the Wigner function [j| . The Husimi 
Q-function can be directly measured by an eight-port ho- 
modyne apparatus or by projection on the harmonic os- 
cillator ground state [10|. These phase-space distribu- 
tions are created to be the closest quantum analogs to a 
classical probability distribution. In this sense, they are 
inherently amenable to direct measurement. 

Weak measurement. — We begin by considering what 
happens to our method for directly measuring the wave- 
function when the state is not pure. At the heart the 
direct method is weak measurement [111 ]. Over the last 
decade, interest in weak measurement has grown as re- 
searchers have realized its potential for interrogating 
quantum systems in a coherent manner [l2| . It has been 
used to model and understand photonic phenomena in 
birefringent photonic crystals [l3j], fiber networks flij . 
cavity QED [15j], and quantum tunneling [lj| 17 1. Weak 
measurement provides insight into a number of funda- 
mental quantum effects, including the role of the uncer- 
tainty principle in the double-slit experiment [3], the 
Legget-Garg inequality [19L the quantum box problem 
[20l | . and Hardy's paradox [2l|. It has also been used to 
amplify small experimental effects 1221 and as feedback 
for control of a quantum system [23j . Weak measure- 
ments have been demonstrated in both classical [13] and 
non-classical systems [25| . 

The concept of weak measurement is universally ap- 
plicable to all types of measurement but here we 
introduce it with a standard model of measurement, the 
von Neumann model [29| . In it, a measurement appa- 
ratus has a pointer, in an initial position wavefunction 
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(sl^i) = ex P [ — I 2 1 (4c 2 )]: whose momentum K 

is coupled with strength g to an system observable A via 
the interaction, U = exp (~igA\£.t/K), where t is the in- 
teraction duration. In a measurement of A, the position 
Q of the pointer is shifted to indicate the result of the 
measurement, A — a: 4>i(q) — > 4>f{q) — 4>i{q — a). In a 
standard (i.e. 'strong', gt large) measurement, this shift 
is much greater than the width a of the pointer and, 
thus, unambiguously indicates the measurement result. 
It also leaves the system in the associated eigenstate |a), 
thereby radically disturbing it. To perform a weak mea- 
surement one reduces gt such that induced pointer shift is 
less then cr, making the measurement result ambiguous. 
The benefit is that the system disturbance is reduced. 
While a weak measurement on a single system provides 
little information, by repeating it on an arbitrarily large 
ensemble of identical systems one can determine the av- 
erage measurement result with arbitrary precision. We 
call this the 'Weak Average' (A w ) p . Unsurprisingly, it is 
simply equal to the standard quantum expectation value 
(A w ) p = Tr [Ap] = (A s ) pl where the latter indicates 
A is measured strongly. 

A distinguishing feature of weak measurement is that, 
in the limit of zero interaction (gt = 0), the quantum 
state of the system remains unchanged. Subsequent mea- 
surements can now provide additional information about 
that initial quantum state \^f). Consider a subsequent 
strong measurement of observable C that results in out- 
come c (corresponding to eigenstate |c)). The average 
result of the weak measurement of A in the sub-ensemble 
of system s gi ving C = c is called the 'Weak Value' and is 
given by [ll[, 
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Surprisingly, the weak value can be outside the range 
of the eigenvalues of A and can even be complex 30- 
[331 ] . Often consideration is limited to its real part, as 
would be done in standard measurement 
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but the imaginary part also has a physical significance: 
the evolution U not only shifts the average position of 
pointer but also the average momentum of the pointer. 
These two, purely real, shifts are proportional to the 
real and imaginary parts of the weak value, respec- 
tively [13-|33: (A w )^, = (Q) / /gt+i (K) f 2a 2 / ' gth, where 
(L)^ = (0/| L I </>/). This result was generalized to other 
initial pointer wavefunctions [Hj], and discrete pointers 
(e.g. Qubits or Spins) [3l|. The complex nature of the 
weak value is what enables us to directly measure the 
real and imaginary parts of the wavefunction and, as we 
show later, directly measure the Dirac distribution and 
density matrix. 

Direct measurement of the quantum wavefunction. — 
We now review our method for the direct measurement of 
the wavefunction. The concept is general, however here 
we consider the case of a discrete Hilbert space. In this 



space, one is free to choose the basis {\a}} (associated 
with observable A) in which the wavefunction will be 
measured. The method consists of weakly measuring a 
projector in this basis 7r Q = \a) (a\, and post-selecting 
on a particular value 60 of the complementary observable 
B. By 'complementary' we mean that (a|&o) = l/VA 
for all a, where N is the dimension of the Hilbert space. 
That is, the overlap is real and constant as function of o. 
The existence of state \bo) is guaranteed by the existence 
of at least two mutually unbiased bases (MUB) in any 
Hilbert space 35j. As discussed in the supplementary 
information of Ref. P, Q , the choice of state bo out of 
the basis { | b) } is simply a convention and is equivalent to 
choosing a reference frame for the direct measurement, 
thereby setting the phases of the basis states in {\a)}. 
Using Eq. ([T]), the quantum state \^>) is given by: |^) = 



v ■ J2 a (^a )* l a )j where (it™ is the weak value and v 
is a constant that is independent of a. Thus, by stepping 
through the values of a in a series of weak measurements 
one can directly measure \^} represented in the a basis. 

Weak measurement of mixed states. — The weak value 
of a system described by a density operator p was first 
considered in [l5j | and shown to be: 



(c|Ap| c ) 
(c\p\c) 



(2) 



Applying this to our direct measurement method we find, 
«) p ° = (bo\a) (a\p\b )/(b \p\b ). Only 2N real pa- 
rameters are found by scanning a. This will not generally 
be sufficient to determine all the parameters in p, which 
has N 2 — 1 real parameters. Consequently, our method 
for the direct measurement of the wavefunction cannot 
be used to determine a mixed state. 

Direct measurement of the Dirac distribution. — We 
now consider what happens if one replaces the strong 
measurement of B with a weak measurement. Specifi- 
cally, we investigate the weak measurement of the prod- 
uct of projectors from the two MUB, S a b = |^)(H a )( a | 
with no post-selection whatsoever. We wish to mea- 
sure its weak average. Although, a non-Hermitian oper- 
ator A is typically considered to be unobservable, later, 
we shall outline specific methods to weakly measure it. 
Surprisingly, we show that even if A is non-Hermitian 
(A w ) p = Tr [Ap] still holds for the weak average. In this 
case, (A w ) p is complex with a physical significance sim- 
ilar to that of the weak value (i.e. shifts in the position 
and momentum of the pointer). For now, we use this 
result to find the weak average of S a &: 



(S: b ) p = Tr [S ab p] = (a\ p \b) (b\a) = S p (a, 6), (3) 

where S p (a,b) is the discrete Hilbert space version of 
the Dirac distribution as defined in 36]. Dirac intro- 
duced this phase-space distribution as a wa y to represent 
a quantum operator O in his 1945 paper [37|, "On the 
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Analogy Between Classical and Quantum Mechanics." In 
various guises it has been investigated periodically dur- 
ing last half-century [3^]. In optics, variations of the 
Dirac distribution have been used widely, appearing in 
Walther's definition of the radiance function in radiome- 
trv l3p and Wolf's specific intensity [!o[ (as pointed out 
in [36[). If O = p, the Dirac distribution is a represen- 
tation of the quantum state of a system. For instance, 
the joint weak measurement of a position x and a mo- 
mentum p (i.e. S xp = \p)(p\x)(x\) on a mixed state p 
gives the phase-space version of the Dirac distribution, 
S p (x,p), which, although it is complex, shares many of 
the desired features of a quasi-probability distribution 

In our weak measurement, if one scans a and b, so 
as to directly measure the Dirac distribution over all 
values of (a, b), one completely determines the density 
operator. But in order to actually calculate the den- 
sity operator from the Dirac distribution one must know 
(6 1 a) = exp (iQab) /vN. Since it is not generally known 
what are the bases in the MUB set (for any given Hilbert 
space) a general formula for 9 a b is also unknown. How- 
ever, if {\d}} is taken to be the standard basis (i.e. 
J2a=o \ a ) H = I> * ne identity operator) then one MUB, 
which we take to be {]&)}, will always be the Fourier 
basis HHim, \b) = J2a=o 1°) ex P {i^ab/N) /VN. In 
this case, 9 a b — —2nab/N, where a and b are inte- 
gers solely used to enumerate the states such that < 
a, b < N — 1. With these choices for our complementary 
bases the density operator is simply related to the Dirac 
distribution by a Discrete Fourier Transform, p ai a 2 — 

EbJo 1 S P (Ol, b) e ^b(a 1 -a 2 )/iV ) where paia2 = (fli | p | fl2) . 

This explicitly shows the weak average, (S^,) , contains 
the same information as the density operator. Much like 
the Wigner function, the Dirac distribution can used to 
find the expectation value of an observable through a 
simple overlap integral 36| . Unlike the Wigner function, 
it is compatible with Bayes' law and, thus, is consistent 
with a quantum analog of classical determinism 42 1. 

Direct measurement of the density operator. — While 
quasi-probability distributions are informationally equiv- 
alent to the density operator, they are less commonly 
known and used. Motivated by a desire to understand 
the nature of the density operator we now describe how 
to measure it directly in a given basis. Consider the weak 
measurement of the product of three projectors, Tl aia2 — 
7r O3 7r& 7r Ol , where ir ai = \ai) (ai|,7r& = \b ) (6q| , and 



Tt a2 = \ a 2) (02) and bo is chosen so that (a\b ) = 1/\N 
for all a (but is not required to be from the Fourier basis). 
As before, there is no post-selection. The weak average 
is: 



(n- Q2 > ={ ai \p\a 2 )/N = p aia jN. 



(4) 



weakly measuring the corresponding projectors, 7T ai and 
7r a2 , in between which is a third weak measurement of 
7Tft . The proportionality constant iV _1 can be eliminated 
through the normalization of the density matrix so that 
Tr [p] = 1. Keeping bo fixed while scanning a\ and a-i 
allows one to map out the entire density matrix. 

Weak measurement of products of complementary vari- 
ables. — One cause for concern in our two direct measure- 
ment methods is that S a b and II aia2 are not Hermitian, 
which, according to the postulates of quantum mechan- 
ics, means they are not observable 0. Indeed, coupling 
such operators to a pointer via the von Neumann inter- 
action (as in U) leads to an unphysical non-unitary evo- 
lution. This issue can be circumvented by dividing the 
measurement into a sequence of unitary von Neumann 
interactions. Each has a pointer beginning in same ini- 
tial state 4>i(q). We now describe a pair of schemes that 
use this strategy to weakly measure the product of two 
non-commuting observables E and F thereby measuring 
their weak average ((EF) w ) p . In the process, we will 
show that the weak average ((EF) w ) p = Tr [EFp]. And 
later we will show that weakly measuring just two observ- 
ables is sufficient to implement both direct measurement 
methods. 

Scheme 1 : The first scheme follows a commonly used 
strategy for standard (strong) measurements: perform 
independent measurements of two observables and cor- 
relate the results to find the observables' product. With 
von Neumann measurements the total evolution is Ut = 
exp (^EK^i/ft) exp (igiFK.it/ h), where the subscripts 
indicate observables on pointer 1 or 2. If E and F com- 



mute, strong measurements give ((EF) 55 



(Q1Q2 



7- 

The weak measurement version of this strategy was pro- 
posed in [44]]. It was simplified in [3l| by forming the 
composite operator = Q/2er + iK.a/h, which has the 
form of an annihilation operator (i.e. a \4>i) — 0), so 
that the standard weak value (Eq. [1]) has the simple 
form (A w )^ = (2a/gt)(a) f . Following [H, one can 

show that in the limit gi<?2 (t/0-) 2 <C 
Ut induces the pointer shifts ((EF) W | 

\2 



1. 



Thus, each element p ai a 2 °f the density matrix in any 
chosen basis (here {|a)}) can be measured directly by 



the evolution 
p = Tr [EFp] = 

gi<?2 (2<r/f) z (aid g) f ■ This scheme was demonstrated ex- 
perimentally in [21| for products of commuting observ- 
ables. Ref. [34| showed that it is valid even for non- 
commuting observables E and F if they are measured 
sequentially, as in Ut- (This result can be generalized 
to an n-product observable, such as the triple product 
n ai a 2 [3ll. l34|.) Thus, just as with strong measurement, 
by performing independent measurements of each observ- 
able and then evaluating a joint expectation value on the 
pointers one can measure ((EF) w ) p . 

Scheme 2: The second scheme measures F and 
then, conditioned on the result, measures E, thereby 
measuring their product. A von Neumann interac- 
tion couples F to a first pointer, which shifts the 
pointer's position. With a strength proportional to 
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this shift a second von Neumann interaction couples 
E to a second pointer. This conditional sequen- 
tial measurement is described by the total interaction 
U D = exp (— ig^EK^Qii/Ti) exp (— ig^FTiit/H), where 
the rightmost interaction couples to either D = Q or 
K, and the subscripts refer to pointer 1 or 2. For 
D = K our weak measurement of F is a standard von 
Neumann interaction. In the limit of gK92t 2 /<J <C 1, 
the evolution Uk shifts the position of pointer 2 by 
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= 0, 



Q 2 )/ = (gK92t 2 ) Re {Tr [EFp]}. However, (K 2 
leaving us without Im{Tr[EFp]}. The imaginary com- 
ponent can be found by coupling F to the position, 
rather than momentum, of the first pointer. If we set 
D = Q, then in the limit gQg2t 2 a/h <C 1, the evo- 
lution Uq shifts the momentum of the second pointer 
by {Q2) f = (2 5q52 t 2 a 2 /n)lm{Tr[EFp]}. In summary, 
conditional two sequential weak measurements allow us 
to measure the real and imaginary components of expec- 
tation value of a product operator. 

Both of the schemes can be followed by post-selection 
of some other observable and, in that case, they would 
give the weak value of EF, which can be complex. How- 
ever, our schemes show that, even without post-selection, 
the weak average of a non-Hcrmitian EF will be com- 
plex. This may come as a surprise since post-selection is 
often cited as the mechanism for anomalous weak values 

How to substitute a strong measurement for a weak 
measurement of one observable in a product. — We now 
show that the weak measurement of CG is equivalent 
to first weakly measuring G followed by a strong mea- 
surement of C, ((CG) W ) = (C S G W ). This can greatly 
simplify both of our proposed methods. In analogy to 
a standard joint expecation value, by (C S G W ) we mean 
that for each measured outcome of the strong measure- 
ment of C one multiplies the corresponding eigenvalue c 
by its probability P(c) by the weak value (G w )^ to find 
the average, £ c cP(c) (G w )" = (C S G W ). Using Eq. ©, 
it follows that, (C S G W ) = £ c c (c\ Gp |c) = Tr [CGp] = 
((CG) W ). In other words, a joint weak-strong measure- 
ment of G and C, respectively, will have the average re- 
sult, Tr [CGp]. Note that because both the weak value 
and corresponding the pointer expectation values (e.g. 
(Q)/) are normalized by -P(c), in an actual experiment 
the pointer signal will be proportional to (c| Gp |c) di- 
rectly, removing the need to find the weak value. 

Consequently, we can directly measure the full Dirac 
distribution and density matrix by measuring the cor- 
relations between the weak measurement and the sub- 
sequent strong measurement outcome while we scan the 
weak measurement. For the density matrix one would 
weakly measure 7T{, 7r ai using scheme 1 or 2 and then 
strongly measure of 7T a2 . Alternatively, one can instead 
strongly measure A, which projects on the eigenstates 
{■7r a } in parallel. This has the advantage that the final 
measurement no longer must be scanned. Similarly, in 



the experiment in [ij, replacing the slit with a camera 
(with a preceding polarization analyzer) in the momen- 
tum plane is the only modification necessary to directly 
measure the Dirac distribution for the transverse den- 
sity matrix of a photon. Evidently, dividing the joint 
measurement into weak and strong parts simplifies each 
of the two new direct measurements considered in this 
paper (given in Eq. ([3|) and Eq. (j4|)) and makes them 
feasible with existing technology (e.g . weak measurement 
of the product of two operators [21j). 

Conclusion. — In this work, we have shown that by 
weakly measuring pairs or triple products of comple- 
mentary variables of a system it is possible not only to 
directly measure its wavefunction but, also, its density 
operator. To measure over the extent of wavefunction 
we only need to scan the first variable of a complemen- 
tary pair. To determine the density operator through its 
Dirac distribution we must, additionally, scan the second 
variable. To determine the density matrix directly, one 
weakly measures the product of a variable, its comple- 
mentary variable, and the first variable again. Leaving 
the complementary variable fixed at one value while ras- 
tering the values of the other two completely maps the 
density matrix, one element at a time. This procedure 
thus provides the density matrix with an operational def- 
inition, it is the average result of a joint weak measure- 
ment of a variable, then its complementary variable, then 
the original variable. 

These methods also provide alternatives to standard 
state tomography that have three key advantages: One, 
they are simpler in that they only require measurements 
in two of the system's bases. Two, they do not require a 
global reconstruction — states can be determined locally, 
point by point. And three, the amount of state distur- 
bance can be minimized. Thus, in principle, we can char- 
acterize quantum states in situ, for instance, in the mid- 
dle of quantum computation circuits, or during chemical 
reactions, without disturbing the process in which they 
feature. 

We thank Aephraim Steinberg for useful discussions. 
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